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Abstract 

The pole mass of a heavy quark is ambiguous by an amount of order Aqcd- We show 
that the heavy-quark potential, V(r), is similarly ambiguous, but that the total static 
energy, 2M po i e + V(r), is unambiguous when expressed in terms of a short-distance 
mass. This implies that the extraction of a short-distance mass from the quarkonium 
spectrum is free of an ambiguity of order Aqcd, in contrast with the pole mass. 



The pole mass of a heavy quark is known to be well defined at any finite order in 
perturbation theory [0]. However, it is known that the pole mass of a heavy quark is an 
ambiguous concept at large orders in perturbation theory || [|, ||. The perturbative series 
relating the pole mass to a short-distance mass (such as the MS mass) has coefficients 
which grow factorially, leading to an ambiguity of order Aqcd in the pole mass || ||. This 
factorial divergence is related to the existence of an infrared renormalon pole || in the 
Borel-transformed quark self-energy when evaluated on shell [|J. 

The ambiguity in the heavy-quark pole mass can be understood heuristically in terms of 
confinement. The pole mass of a heavy meson (bound state of a heavy quark and a lightf] 
antiquark) is well defined, since it is a physical quantity. To extract the pole mass of the 
heavy quark from the meson mass, one must subtract the binding energy, of order Aqcd- 
However, it is impossible to define this binding energy unambiguously, since one cannot 
separate the light antiquark from the heavy quark, due to confinement. Hence the pole mass 
of a heavy quark is ambiguous by an amount of order Aq C d- 

Consider instead quarkonium, a bound state of a heavy quark and a heavy antiquark. If 
the mass of the heavy quark is sufficiently large, the strong coupling governing the heavy- 
quark interactions in quarkonia states with small principal quantum number is weak ||. 
The quarkonium can then be handled just like positronium in electrodynamics, by sum- 
ming uncrossed Coulomb ladder diagrams to yield the Schrodinger equation, and using non- 
relativistic perturbation theory. Thus the low-lying quarkonia states can be described almost 
entirely in terms of QCD perturbation theory. It is tempting to conclude that, if the quark 
mass is heavy enough, the pole mass can be extracted to arbitrary precision from a per- 
turbative calculation of the binding energy of these states, perhaps supplemented by some 
nonperturbative input. The goal of this article is to show that this is not the case; the heavy- 
quark pole mass cannot be extracted from the quarkonium spectrum with an accuracy better 
than order Aqcd- This is an important point to clarify, since the c- and 6-quark pole masses 
are extracted from the J/ip 0, |, g, |I| [n] and T spectra Q, |, |, [jj, [Ej| |16| . 



We begin by considering the static potential between a heavy quark and a heavy antiquark 
in a color singlet state at leading order in QCD |17|, [18], [19], |2(| 



J^-^ = - Cf T (^ = 4 /3) , (1) 

where a s = a s (n) is the strong coupling evaluated at the scale /i. An elegant means to 
analyze large orders in perturbation theory is to calculate the Borel transform (with respect 
to boa s / A.n) of the static potential. This calculation, performed in Ref. p2| , givesQ 

d 3 k e ikr 



V(r) 




(27r) 3 k 2 ( 1 +«) 



b ' r(2u + 1) 

C^e-^^-^-r- + • • • , (2) 
bo [u - 2) 



1 m < Aqcd 

2 An analysis of quarkonia energies at large orders in perturbation theory based on the position-space 
potential V(r) — —Cf as ^ r ^ may be found in Ref. |p3[. 
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where u is the Borel parameter, C is a renormalization-scheme-dependent constant (C = 
—5/3 in the MS scheme), and 

&o = ll-~iV/ (3) 

is the one-loop beta- function coefficient. The infrared renormalon pole nearest the origin, 
at u = 1/2, controls the asymptotic behavior of the perturbation series. This pole and its 
residue are made explicit in the last line of Eq. (0) . 

The static potential is recovered from Eq. @ by inverse Borel transformation, 

poo _ 

V(r)= due~ Anu/{baas) V{r) . (4) 



The evaluation of this integral is impeded by the presence of infrared renormalon poles in 
V(r) on the positive real axis at all half-integer values of u. The ambiguity in the integral is 
dominated by the infrared renormalon pole closest to the origin, at u = 1/2. Estimating the 
ambiguity as half the magnitude of the difference between deforming the integration contour 
above and below the pole yields 

SV (r) ~ CV-^e-^V- 2 ^" ^ 
bo 

~ C F %^A QCD , (5) 
bo 

where the one-loop renormalization-group equation for a s has been used to obtain the final 
expression. The ambiguity is a constant shift of the potential, by an amount of order Aqcd- 
The ambiguity of order Aqcd in the static potential was first derived in Ref. |22]| , but no 
attempt was made to interpret it. We now show that this ambiguity is due to the (implicit) 
use of the pole mass, which is sensitive to momenta of order Aqcd, in the definition of the 
static potential. 

The relation between the heavy-quark pole mass, M po i e , and a short-distance mass, M, 
is given in Borel space by || 

1 o \ /i / 1 (3 — u) u 




= -CfU- CI ^^^ + ■ ■ ■ , (6) 
bo 

where R-e 1 (u) is a renormalization-scheme-dependent regular function. The infrared renor- 
malon pole closest to the origin is again at u = 1/2, and its residue is displayed in the last 
line of Eq. (H). The ambiguity in the pole mass can be estimated as above for the static 
potential, and is given by 0] 

6M pole ~ C F ^e- c/2 A QC D , (7) 
bo 

as is well known. 
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Consider now the total static energy of a heavy quark and a heavy antiquark in a color 
singlet state with spatial separation r, given by the sum of the static potential energy and 
the rest mass of the particles. In Borel space, 

E sta tic(r) = 2M pole + V(r) . (8) 

If the pole mass is eliminated in favor of a short-distance mass via Eq. we see that 
the infrared renormalon pole at u = 1/2 in the total static energy is cancelled. Thus the 
total static energy does not have an ambiguity of order Aqcd when expressed in terms of a 
short- distance mass. 

The total static energy of the heavy quark-antiquark pair in a color-singlet state is an 
unambiguous concept ^| [T7|]. However, when one attempts to separate this energy into the 
sum of the static potential energy and the pole masses of the quark and antiquark |L7| [18], [19], 
|2"0| , |2T[| , the ambiguity of order Aqcd in the pole mass results in a corresponding ambiguity in 
the static potential. This is the source of the ambiguity of order Aqcd in the static potential. 
The ambiguity disappears when the pole mass is replaced with a short-distance mass. This 
well-known phenomenon occurs for many other processes in QCD [|[ |j, ^5], ^6], ^7], ^9], pO 



Because it will be useful to us when we consider the full dynamical (nonstatic) quarkonium 
calculation, let us recall that an analysis of the infrared behavior can also be carried out via 
a finite gluon mass A || |25|, |3l| . The static potential becomes 



V(r) = -C F —e~ Xr = -C F a s ( - - A + ■ • ■ ) . (9) 

r \r / 

The term of order A corresponds to the infrared renormalon pole at u = 1/2. The relation 



between the pole mass and a short- distance mass is P, 25, ET\ 



M pole = M-C F ^-X + --- . (10) 

We see that the term linear in the gluon mass cancels when we calculate the total static 
energy, 2M po i e + V(r), in terms of a short-distance mass. This corresponds to the cancellation 
of the pole at u = 1/2 in the Borel transform of the total static energy, Eq. (|8]). 

The full dynamical (nonstatic) quarkonium calculation at leading order in the nonrela- 
tivistic approximation requires solving the Schrodinger equation 



V 
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+ V(r)-E)G(r,0,E) = 5^(r), (11) 

where G(r, 0, E) is the Schrodinger-equation Green function, and E = — 2M po i e is the 
binding energy. The total center-of-mass energy y 7 ? is physical and unambiguous, but the 
binding energy E is not, since it requires subtracting twice the pole mass. Eliminating E in 
favor of yfs in Eq. ([□]) gives 

/ V ~ + 2M poie + l/(r) - ^ ] G(r, 0, ^ - 2M poie ) = 5( 3 »(r) . (12) 



M po i e 



The quarkonia masses correspond to the poles of the Green function in the ^/s plane. We 
see that the total static energy, 2M po i e + V(r), now appears in the Schrodinger equation. As 
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shown above, the total static energy is not ambiguous by an amount of order Aqcd, but due 
to the ambiguity in V(r), the pole mass cannot be extracted from the quarkonium spectrum 
with an accuracy better than order Aqcd- However, the static potential can be made free of 
the ambiguity of order Aqcd if the pole mass is exchanged for a short- distance mass and if 
twice the difference between the pole and the short-distance mass is absorbed into the static 
potential, 

f + 2M + V(r) - ] G(r, 0, ^ _ 2 M) = 5^(r) , (13) 

\ Mpoie J 

where 

V{r) = V(r) + 2(M pole - M) . (14) 

Hence the accuracy with which a short-distance mass can be extracted from the quarkonium 
spectrum is not limited by order Aq C d- 

The pole mass also appears in the denominator of the kinetic-energy term in the Schro- 
dinger equation. However, replacing this mass with a short-distance mass only affects terms 
suppressed by powers of a s . This can be seen most easily by using the gluon mass as an 
infrared regulator, Eqs. (P) and (|H]). The linear gluon mass terms which cancel in the total 
static energy are of order a s X. Since the kinetic-energy term is of order M po i e a 2 s , the linear 
gluon mass term generated by replacing M po i e with a short- distance mass is of order a^X. 

In this paper we have shown that the heavy-quark pole mass cannot be extracted from 
the quarkonium spectrum with an accuracy better than order Aqcd- This is relevant for 
the determination of the c- and 6-quark pole masses from the J /if) J7|, [| ||, [10], |TTJ and T 
spectra j7|, §, || [12|, [13], [14], [15], |16| . However, the accuracy with which a properly-defined 



short-distance mass can be extracted from the quarkonium spectrum is not limited by order 
Aqcd- 

Note added: The results of this paper have also been arrived at in a recent paper 



Acknowledgments 

We are grateful for conversations with A. El-Khadra. A. H. H. was supported in part 
by the U. S. Department of Energy under contract No. DOE DE-FG03-90ER40546. M. S., 
T. S., and S. W. were supported in part by the U. S. Department of Energy under contract 
No. DOE DE-FG02-91ER40677. M. S. was supported in part by a grant from the UIUC 
Campus Research Board. 

References 

[1] A. Kronfeld, Phys. Rev. D 58, 051501 (1998). 

[2] I. I. Bigi and N. G. Uraltsev, Phys. Lett. B321, 412 (1994). 

[3] I. I. Bigi, M. A. Shifman, N. G. Uraltsev, and A. I. Vainshtein, Phys. Rev. D 50, 2234 
(1994). 



4 



[4] M. Beneke and V. M. Braun, Nucl. Phys. B426, 301 (1994); 
M. Beneke, Phys. Lett. B344, 341 (1995). 

[5] G. 't Hooft, in The Whys of Suhnuclear Physics, Proceedings of the International School 
of Subnuclear Physics, Erice, 1977, ed. A. Zichichi (Plenum, New York, 1979), p. 943; 
A. Mueller, in QCD - 20 Years Later, Proceedings of the Workshop, Aachen, Germany, 
1992, eds. P. Zerwas and H. Kastrup (World Scientific, Singapore, 1993), Vol. 1, p. 162. 

[6] T. Appelquist and H. D. Politzer, Phys. Rev. Lett. 34, 43 (1975). 

[7] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 

[8] S. Narison, Phys. Lett. B197, 405 (1987); Phys. Lett. B341, 73 (1994); 
Phys. Lett. B352, 122 (1995). 

[9] S. Titard and F. J. Yndurain, Phys. Rev. D 49, 6007 (1994); 
A. Pineda and F. J. Yndurain, Phys. Rev. D 58, 094022 (1998). 

[10] C. Dominguez, G. Gluckman, and N. Paver, Phys. Lett. B333, 184 (1994). 

[11] A. Kronfeld, Nucl. Phys. Proc. Suppl. 63, 311 (1998). 

[12] C. Davies et al, Phys. Rev. Lett. 73, 2654 (1994). 

[13] M. B. Voloshin, Int. J. Mod. Phys. A10, 2865 (1995). 

[14] M. Jamin and A. Pich, Nucl. Phys. B507, 334 (1997). 

[15] J. H. Kiihn, A. A. Penin, and A. A. Pivovarov, Nucl. Phys. B534, 356 (1998); 
A. A. Penin and A. A. Pivovarov, Phys. Lett. B435, 413 (1998). 

[16] A. H. Hoang, Phys. Rev. D 59, 014039 (1999). 

[17] L. Susskind, in Weak and Electromagnetic Interactions at High Energy, Les Houches 
1976, eds. R. Balian and C. Llewellyn Smith (North Holland, Amsterdam, 1977), p. 207. 

[18] W. Fischler, Nucl. Phys. B129, 157 (1977). 

[19] T. Appelquist, M. Dine, and I. Muzinich, Phys. Lett. B69, 231 (1977). 
[20] F. Feinberg, Phys. Rev. Lett. 39, 316 (1977). 

[21] M. Peter, Phys. Rev. Lett. 78, 602 (1997); Nucl. Phys. B501, 471 (1997). 
[22] U. Aglietti and Z. Ligeti, Phys. Lett. B364, 75 (1995). 
[23] M. C. Smith, T. Stelzer, and S. Willenbrock, Phys. Lett. B413, 405 (1997). 
[24] K. Wilson, Phys. Rev. D 10, 2445 (1974). 

[25] M. Beneke, V. M. Braun, and V. I. Zakharov, Phys. Rev. Lett. 73, 3058 (1994). 



5 



[26] M. Beneke and V. M. Braun, Phys. Lett. B348, 513 (1995). 

[27] B. H. Smith and M. B. Voloshin, Phys. Lett. B340, 176 (1994). 

[28] P. Gambino and A. Sirlin, Phys. Lett. B355, 295 (1995). 

[29] M. Luke, A. Manohar, and M. Savage, Phys. Rev. D 51, 4924 (1995). 

[30] M. Neubert and C. Sachrajda, Nucl. Phys. B438, 235 (1995). 

[31] P. Ball, M. Beneke, and V. M. Braun, Nucl. Phys. B452, 563 (1995). 

[32] M. Beneke, Phys. Lett. B434, 115 (1998). 



6 



